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Abstract 

We propose an alternate, calculable mechanism of dark matter genesis, "thermal freeze- 
in," involving a Feebly Interacting Massive Particle (FIMP) interacting so feebly with the 
thermal bath that it never attains thermal equilibrium. As with the conventional "thermal 
freeze-out" production mechanism, the relic abundance reflects a combination of initial ther- 
mal distributions together with particle masses and couplings that can be measured in the 
laboratory or astrophysically. The freeze-in yield is IR dominated by low temperatures near 
the FIMP mass and is independent of unknown UV physics, such as the reheat temperature 
after inflation. Moduli and modulinos of string theory compactifications that receive mass 
from weak-scale supersymmetry breaking provide implementations of the freeze-in mecha- 
nism, as do models that employ Dirac neutrino masses or GUT-scale-suppressed interactions. 
Experimental signals of freeze-in and FIMPs can be spectacular, including the production 
of new metastable coloured or charged particles at the LHC as well as the alteration of big 
bang nucleosynthesis. 
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1 Introduction 

Many theories of Dark Matter (DM) genesis are based upon the mechanism of "thermal freeze out" 
[1] . In this process DM particles have a large initial thermal density which, as the temperature of 
the hot plasma of the early universe drops below the mass, dilutes away until the annihilation to 
lighter species becomes slower than the expansion rate of the universe and the comoving number 
density of DM particles becomes fixed. The larger this annihilation cross section the more the DM 
particles are able to annihilate and hence a thermal distribution with an exponential Boltzmann 
factor is maintained to a lower temperature, giving a lower final yield. An attractive feature 
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x = m/T 



Figure 1: Log-Log plot of the evolution of the relic yields for conventional freeze- 
out (solid coloured) and freeze-in via a Yukawa interaction (dashed coloured) as a 
function of x = m/T. The black solid line indicates the yield assuming equilibrium is 
maintained, while the arrows indicate the effect of increasing coupling strength for the 
two processes. Note that the freeze- in yield is dominated by the epoch x ~ 2 — 5, in 
contrast to freeze-out which only departs from equilibrium for x ~ 20 — 30. 

of the freeze-out mechanism is that for renormalisable couplings the yield is dominated by low 
temperatures with freeze-out typically occurring at a temperature a factor of 20 — 25 below the 
DM mass, and so is independent of the uncertain early thermal history of the universe and possible 
new interactions at high scales. 

Are there other possibilities, apart from freeze-out, where a relic abundance reflects a com- 
bination of initial thermal distributions together with particle masses and couplings that can be 
measured in the laboratory or astrophysically? In particular we seek cases, like the most attractive 
form of freeze-out, where production is IR dominated by low temperatures of order the DM mass, 
m, and is independent of unknown UV quantities, such as the reheat temperature after inflation. 

In this paper we show that there is an alternate mechanism, "freeze-in", with these features. 
Suppose that at temperature T there is a set of bath particles that are in thermal equilibrium and 
some other long-lived particle X, having interactions with the bath that are so feeble that X is 
thermally decoupled from the plasma. We make the crucial assumption that the earlier history 
of the universe makes the abundance of X negligibly small, whether by inflation or some other 
mechanism. Although feeble, the interactions with the bath do lead to some X production and, 
for renormalisable interactions, the dominant production of X occurs as T drops below the mass 
of X (providing X is heavier than the bath particles with which it interacts). The abundance of 
X "freezes-in" with a yield that increases with the interaction strength of X with the bath. 

Freeze-in can be viewed as the opposite process to freeze-out. As the temperature drops below 
the mass of the relevant particle, the DM is either heading away from (freeze-out) or towards 
(freeze-in) thermal equilibrium. Freeze-out begins with a full T 3 thermal number density of DM 
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particles, and reducing the interaction strength helps to maintain this large abundance. Freeze- 
in has a negligible initial DM abundance, but increasing the interaction strength increases the 
production from the thermal bath. These trends are illustrated in Figure [TJ which shows the 
evolution with temperature of the dark matter abundance according to, respectively, conventional 
freeze-out, and the freeze-in mechanism we study here. 

In Section |2j as well as outlining the basic freeze-in mechanism and comparing its features 
with those of freeze-out, we also introduce the idea of a FIMP — a "Feebly- Interacting-Massive- 
Particle" (or alternatively "Frozen-In-Massive-Particle" ) — as distinct from a WIMP whose relic 
abundance is set by conventional freeze-out. Two cases are considered: either the FIMP itself is 
the dark matter which is frozen-in, or the dominant contribution to the DM density arises from 
frozen-in FIMPs which then decay to a lighter DM particle. For enhanced pedagogy the detailed 
calculation of the DM abundance in these cases is postponed until Section [6] 

We turn in Section [3] to the question of motivated DM candidates that have a relic abundance 
determined by the IR-dominated freeze-in mechanism, and show that the moduli and modulinos 
of compactified string theories with weak-scale supersymmetry breaking provide implementations 
of the freeze-in mechanism, as does any extra-dimensional extension of the Standard Model (SM) 
with some moduli stabilised at the weak scale. Models following from Dirac neutrino masses, 
or involving kinetic mixing with a hidden sector also naturally accommodate FIMPs and the 
freeze-in mechanism. A striking difference with freeze-out production of DM, is that for freeze-in 
the final relic density is, in the simplest cases, automatically independent of the FIMP mass, 
allowing superheavy as well as weak-scale mass DM candidates, and we touch on a high-scale 
extra-dimensional realisation of such a scenario. Theories utilising GUT-scale-suppressed non- 
renormalisable interactions involving SM or MSSM Higgs fields which become renormalisable when 
the Higgs gets its vacuum expectation value also naturally give rise to FIMPs. However, in this 
case, there is also a UV contribution to the freeze-in yield which can dominate the IR contribution 
if the reheat temperature is sufficiently large, and we therefore postpone the discussion of such 
implementations until Section [7] where the effect of higher-dimension-operators on freeze-in yields 
is considered. 

Experimental signals of freeze-in and FIMPs are briefly summarised in Section |4j with a com- 
panion paper [2] more extensively discussing the rich phenomenology that accompanies freeze-in 
to follow. The signals depend on the nature of the Lightest Observable Sector Particle (LOSP) 
that carries the conserved quantum number that stabilises DM. (In theories with weak scale su- 
persymmetry, the LOSP is the lightest superpartner in thermal equilibrium at the weak scale.) 
There are two general possibilities: The DM particle is the FIMP, and there are spectacular LHC 
signatures arising from the production of coloured or charged LOSPs, which stop in the detector 
and decay with a long lifetime. In addition, LOSPs decaying to FIMPs in the early Universe 
may lead to interesting modifications of Big-Bang Nucleosynthesis (BBN). Alternatively, the DM 
particle is the LOSP but its interactions are too strong to have a sufficiently large freeze-out relic 
density; instead its relic abundance results from freeze-in of FIMPs, with the FIMPs later decaying 
to the DM. The second scenario can give signals via alteration of big-bang nucleosynthesis element 
abundances, and via increased DM pair annihilation relevant for indirect detection of DM. 

Thermal relic abundances conventionally arise by decoupling of a species that was previously 
in thermal equilibrium, whether with or without a chemical potential. Freeze-in provides the 
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only possible alternative thermal production mechanism that is dominated by IR processes. In 
Section [5] we sketch "abundance phase diagrams" that show regions of mass and coupling where 
each mechanism dominates the production of the relic abundance. The topology of these phase 
diagrams, as well as the number of domains where different mechanisms dominate, depends on 
the form of the DM-bath interaction, as we illustrate with two examples arising from a quartic 
scalar interaction and a Yukawa interaction. Furthermore, we argue that there exists a certain 
"universality" to the phase diagram behaviour at small coupling. 

After presenting the detailed calculation of the relic abundance in various cases in Section |6j 
we discuss the physics of FIMP interactions mediated by higher-dimension-operators, as well as 
some variations of the basic freeze-in mechanism in Section [7| Finally we conclude in Section [8j 



2 General Features of Freeze-In 

The basic mechanism of freeze-in is simple to describe although, as we will argue later, there can 
be many variations with important differences of detail and signals. Here we give the general 
mechanism. At temperatures well above the weak scale we assume that there is a FIMP, X, 
that is only very weakly coupled to the thermal bath via some renormalisable interaction. The 
interaction vertex may involve more than one particle from the thermal bath and the mass of the 
heaviest particle at the vertex is m, which we typically take to be near the weak scale. For a 
Yukawa or quartic interaction, the dimensionless coupling strength is A, while for a trilinear scalar 
interaction the coupling is Am, and in all cases A< 1. 

At very high temperatures we assume a negligible initial X abundance. As the universe evolves 
X particles are produced from collisions or decays of bath particles, but at a rate that is always 
suppressed by A 2 . During a Hubble doubling time at era T ^> m, the X yield is 

2 Mpi ( ^ m 2 \ 



Y{T) ~ X'^f \1, — 2 y (2.1) 

where Y = n/S, n is the number density of X, and S is the entropy density of the plasma. The 
m independent yield corresponds to quartic interactions, while the additional m 2 /T 2 suppression 
applies to Yukawa interactions and the super-renormalisable case. The process is always IR 
dominated, favouring low temperatures. The dominant production occurs at T ~ m, since at lower 
temperatures there will be an exponential suppression resulting from the necessity of involving a 
particle of mass m > T. Hence, for all renormalisable interactions, the abundance of X "freezes-in" 
with a yield 

Y FJ ~ A 2 (MH) . (2.2) 

As mentioned in the Introduction, freeze-in can be viewed as the opposite process to freeze-out. 
We recall that, in the absence of a chemical potential, the freeze-out yield is given by 

Y fo ^ - r (2.3) 
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Figure 2: Schematic picture of the relic abundances due to freeze-in and freeze-out as 
a function of coupling strength. The way in which the freeze-out and freeze-in yield 
behaviours connect to one another is model-dependent. As we show in detail in Section 
[5j freeze-in and freeze-out are in fact two of the four basic mechanisms for thermal DM 
production, and we sketch the "abundance phase diagrams" of DM depending upon 
the strength and type of the DM-thermal bath interaction and the DM mass. 



In the simple case that the DM mass ml is the only mass scale entering av, we have av ~ X' 2 /m' 2 , 
where A' is the relevant interaction strength, giving a freeze-out yield of 

y ™ ~ M£> • {2A) 

Freeze-in assumes a negligible initial X abundance, but increasing the interaction strength, A, 
increases the production from the thermal bath, while freeze-out begins with a full T 3 thermal 
number density of DM particles, and reducing the interaction strength, A', helps to maintain this 
large abundance (see Figure [2]). Indeed, the yields eqns.(2.2) and (2.4) show inverse dependences 
on the coupling and mass, which is stressed by writing 

Y FI ~ X 2 mt m , Y FO ~ - 2 — , (2.5) 

X u m't m i 

where t m ~ Mpi/m 2 (t' m ~ Mpi/m' 2 ) is the Hubble time at the epoch of freeze-in (freeze-out). The 
freeze-out abundance decreases with t m > while the freeze-in abundance is increased by occurring 
at late times. 

Despite these opposite features, freeze-out and freeze-in share crucial common aspects: the final 
out-of-equilibrium abundance, given the relevant particle masses and couplings, can be computed 
solely from an initial state of bath particles that are in thermal equilibrium, and the resulting 
abundance is dominated by IR physics. 

For freeze-out the special case A' ~ 1 and ml ~ v, the scale of weak interactions, gives DM as 
"Weakly Interacting Massive Particles", or WIMPs, with 

Yfo ~ (2-6) 
Mpi 

In practice the cross section may involve more than one mass scale in the TeV domain, so that there 
are orders of magnitude spread in the abundance expected from WIMP dark matter. Nevertheless, 
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the parametric understanding of the abundance of DM is suggestive, and the prediction of a TeV 
mass particle with coupling strengths of order unity offers the hope of collider verification of the 
production mechanism. For freeze-in, if the particle masses are again at the weak scale, the same 
dependence of the relic abundance on v and Mp\ follows if the small coupling A is linear in the 
weak scale: 

* ~ T7- Y FI ~ T7 — ■ (2.7) 

M Pl F/ M Pl v ; 

As in the WIMP case, this parametric behaviour is significantly modified by numerical factors; 
nevertheless, it suggests seeking theories where small couplings arise at linear order in the weak 
scale. 

Whether produced by freeze-out or freeze-in, stable DM is the lightest particle transforming 
non-trivially under some unbroken symmetry. For conventional freeze-out, this lightest particle is 
automatically the WIMP, whereas for freeze-in two particles are of interest: the FIMP and the 
lightest particle in the thermal bath that carries the symmetry, the LOSP. If the LOSP is lighter, 
then LOSP DM is produced by FIMP decay. If FIMPs are lighter, then collider signals involve 
the production of LOSPs followed by decays to FIMPs. In either case, the freeze-in mechanism 
always introduces particles of very long lifetime 

^fimp/losp ~ A m ~ (2.8) 



pi 



where m is the mass of the decaying FIMP or LOSP, and in the last expression we took m ~ v 
and A ~ v/Mpi. The heavier of the LOSP and FIMP will generically decay to the lighter with a 
lifetime within a few orders of magnitude of a second. 

We show in detail in Section [6] that the freeze-in density is dominated, where possible, by 
decays or inverse decays involving the bath particles and X. Freeze-in of a stable X via the 



decays of LOSPs in the thermal bath (see eq. (6.10) and surrounding discussion) gives 



W = 1My ^ mx l B \ (2.9) 



where Bi is the LOSP. For a LOSP decay rate given by Yb x = \ 2 mB 1 /87r and, for simplicity, 
making the good approximation gf ~ g%, the required DM density Vlxh 2 ~ 0.106 occurs for a 
coupling of size 



As discussed in Section |6| this value is reduced to A ~ 1.5 x 10~ 13 (10 2 /gb at / l ) 1/ ' 2 if the FIMP 
couples with comparable strength to many (gbath) bath particles of comparable mass, as occurs, 
for example, in supersymmetric realizations of FIMPs as discussed in Section [3] Thus we learn 
the interesting fact that A ~ v/Mqut works well as an explanation of frozen-in DM. The origin of 
the 1/Mgut suppression could not only be a true "Grand Unified Theory" scale in the traditional 
SU(5) or 50(10) sense, but could alternatively be the string compactification scale Mkk (or 
orbifold GUT compactification scale) or the string scale M str depending on the precise model in 
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which FIMPs are implemented. Note that the choice of M KK ~ (10 15 - 10 16 ) GeV and M str ~ 10 16 
GeV is well motivated by modern string model building [3] pi] [5], or by orbifold GUTs [6]. 

In the following section we argue that FIMP candidates with weak-scale masses and couplings 
of this order naturally appear in many beyond-the-standard-model theories, such as string the- 
ory with low-energy supersymmetry. Moreover in Section [4] we show that such FIMP theories 
can be tested and explored by experiments at the LHC and by astrophysical and cosmological 
observations. 



3 FIMP Candidates 

We now turn to the issue of possible FIMP candidates in well-motivated theories of beyond-the- 
Standard-Model physics. 

3.1 Moduli with weak scale supersymmetry 

The FIMP may be a modulus, T, or modulino, T, frozen-in by a variety of possible supersymmetric 
interactions at the weak scale. We consider a theory with weak scale supersymmetry where the size 
of the supersymmetry breaking arises from the compactification of extra spatial dimensions, and 
depends on T. We assume that the conventional moduli problem is solved, so that after reheating 
at the end of inflation the T abundance is negligible. The modulus T has a large vev, which is 
not determined until supersymmetry is broken, so that T is R parity even and the corresponding 
modulino, T, is R parity odd. The modulus and modulino are expected to pick up masses of order 
the supersymmetry breaking scale, which we take to be of order the weak scale v. It is therefore 
natural to imagine that the scalar T can play the role of the FIMP that is produced by freeze-in 
and then decays to LOSP DM, while the .R-parity odd modulino T can either similarly be freeze-in 
produced and then decay to LOSP DM, or be FIMP DM itself, depending upon the hierarchy of 
LOSP and modulino masses. 

Most important, the interactions of T and T are naturally of the right size to lead to freeze- in 
production. These interactions are obtained by expanding the supersymmetry breaking parame- 
ters about the modulus vev. Taking these parameters to be soft scalar masses, m 2 , the conven- 
tionally defined A and B parameters, the gaugino masses m§ and the /i parameter, the leading 
(renormalisable) interactions of the modulus are 

m2 + If) (0t0 + k]k) ^ i 1 + If) h2 Ay i 1 + If) 02/1 

m- 9 (l + ^) gg ,y (l + J) A* , (l + £) hh, 

where now T refers to the fluctuation of the modulus about its expectation value. Here is a 
squark or slepton, h a Higgsino, and h a Higgs boson, and y is a corresponding Yukawa coupling. 
The mass scale M is the compactification scale, which we take to be of order the unified mass 
scale, M u , as in many string constructions, and numerical factors of order unity are ignored. For 
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the modulino the leading interactions are 



T ~ 

Vjjhh, (3.1) 

as well as possible terms that arise from higher-dimension D-terms involving the moduli-dependent 
susy-breaking spurion leading to interactions with sfermion-fermion pairs 

m 6 



"susy 



f(qq\ll\hh ] ). (3.2) 



These interactions are all of the form shown in eq. (5.4), with 

A ~ ± (3.3) 

in each case. Thus we discover that the modulus and modulino corresponding to the size of 
supersymmetry breaking automatically lead to FIMP-like behaviour. This is not the case for 
moduli that appear in the expansion of renormalisable couplings, since then the moduli couplings 
are non-renormalisable and do not lead to IR-dominated freeze-in behaviour, as we discuss in 
Section [7j 

There are many possible reactions for the freeze-in mechanism involving T or T, and here we 
give a few cases for illustration. T FIMP dark matter could arise from the freeze-in mechanism 
by decays of a LOSP chargino 

X + ^H+f, (3.4) 

where H + is a charged Higgs boson. On the other hand, if the T FIMP is heavier than the LOSP, 
then DM may arise by freeze-in of T followed by the decay of T to LOSP dark matter, for example 

f-tf, f^hx\ (3.5) 

where t and t are the top quark and top squark, and x° the neutralino LOSP. Alternatively, 
neutralino DM could be produced by the freeze-in of a modulus FIMP which then decays, for 
example 

ia-txT, (3-6) 



3.2 Dirac neutrino masses within weak scale supersymmetry 

As a second example with weak-scale supersymmetry, consider the case that neutrinos are Dirac, 
with masses generated via the interaction 

XLNH U , (3.7) 

where L, N and H u are chiral superfields containing the lepton doublet, right-handed neutrino 
singlet and Higgs fields. The coupling matrix A now determines the neutrino masses and conse- 
quently is constrained by experiment to have entries with maximal values that are of order 10~ 13 , 
very close to that required for a FIMP. The interaction leads to the freeze-in of the right-handed 
sneutrino, vr, for example via hi — > ur an d, as computed by Asaka et al [7J, can lead to a successful 
dark matter abundance when vr is the LSFQ 

1 For related ideas see [8] and for the sterile neutrino case see e.g. [9]. 
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3.3 FIMPs from kinetic mixing 



Additional hidden U(l) factors under which no light SM states are charged are a common feature 
of theories beyond-the-Standard-Model. For example, many string constructions possess multiple 
hidden U(l) factors arising from D-branes or from antisymmetric tensor fields with multiplicity 
determined by the topology of the compactification |10j . The unique renormalisable coupling of 
such hidden U (l)'s to SM states arises via the kinetic mixing term, e iY Fj J ' 1 ' F Y ^ U , with hypercharge, 
and in field theory such kinetic mixing can be generated at the one-loop level when super-massive 
states charged under both £7(1) 's are integrated out [UJ, with result 

e iY * ^yQrgiQi) log ( m\ (3 8) 



127T 2 

For small mass splittings between the superheavy states e^y ~ Na/(37r)(Am/m) where N counts 
the multiplicity of heavy states, Am is their mass splitting, and m their mean mass. Although 
these mixings can be often harmlessly rotated away for the photons (as the string £7(1) 's often 
possess no light charged states), they imply non-trivial couplings to the MSSM states of the massive 
hidden photini superpartners via mixing with the bino, leading to hidden photini as candidate 
FIMPs [121 M- 

A natural possibility is that the massive states sit at the GUT, or compactification scale with 
splittings arising from electroweak-symmetry-breaking effects. In this case the field theory calcu- 



lation eq. (3.8) gives e^y ~ 1CT 12 — 1CT 15 , with the higher values applying if the compactification, 
or effective GUT scale, occurs at 10 14 — 10 15 GeV as for orbifold GUTs and many semi-realistic 
string compactifications. For both the heterotic string [13] and type-II string theory [15] kinetic 
mixing with U (l)y arises by a process that generalises this field theory calculation. For example in 
the type-II case stretched open string states with one end on the SM brane stack and the other on 
the hidden D-brane lead to massive states charged under both [/(l)'s, and a one loop open string 
diagram then, in general, generates kinetic mixing [15J. The resulting mixing is model dependent, 
but is exponentially suppressed values if the compactification is warped, or if the mediating fields 
are massive, e.g., due to fluxes. For either case if e iY ~ 1CT 12 — 1CT 15 the hidden photini can be 
FIMPs. 



3.4 Very heavy FIMPs and extra dimensions 



Dark matter with relic abundance Y and mass m leads to a temperature for matter-radiation 
equality that has the parametric form T e ~ Ym. Hence, from eqns.(2.4) and (2.2) freeze-out and 
freeze-in lead to 



T, 



Tfl 

e,FO ~ x /9 „ J , Te,FI ~ A M Pi . (3.9) 

X' z Mpi 

For freeze-out, as m' increases A' must also increase to maintain the observed T e . As is well 
known, unitarity prevents a WIMP being much heavier than the weak scale [TB|. The situation 
is completely different for freeze-in; remarkably T e is independent of the FIMP mass, which can 
therefore take any extremely large value, being bounded only by the reheat temperature after 



r> 
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inflation. Furthermore, the very feeble coupling of the FIMP is a generic feature of freeze-in 

a ~(0 /2 - < 3 - io > 

Suppose there is some stable particle, X, much heavier than the weak scale. If it has a coupling 
to SM particles, for example via a quartic scalar interaction XX^ X H , that allows it to reach 
thermal equilibrium then it overdoses the universe. To avoid this, the coupling must be reduced; 
in fact, it will overdose the universe until the coupling is reduced far enough that it satisfies 
eq. (3.10) and the heavy stable particle becomes FIMP dark matter. How might such a small 
interaction arise? If the SM is localised on a brane in a higher dimensional manifold while the 
heavy particle X lives in the bulk, then a small coupling can arise from the X profile having a 
small value at the location of the SM brane. Very small couplings occur easily since the profiles 
are frequently exponential or Gaussian. 

In the multiverse there may be an environmental requirement that the DM abundance not 
be much larger than we observe. For example, axionic DM with a very large decay constant 
typically overdoses the universe, but this may be avoided by environmental selection of a small 
initial vacuum misalignment angle of the axion field p2], yielding axion DM. Similarly, the bulk 
mass of the very heavy stable particle X may be selected to give a sufficiently steep profile that 
the coupling A obeys eq. (3.10), thus creating FIMP dark matter. 



4 Experimental Signatures of Freeze-In & FIMPs 

No matter what the underlying theory for freeze-in, the coupling of the FIMP to the thermal 
bath is very small, so a crucial question is whether the freeze-in mechanism can be tested by 
measurements at accelerators or by cosmological observations. In this section we outline several 
possible signals although we leave the details to a companion paper [2] . In addition to the FIMP, 
the freeze-in mechanism typically requires a LOSP and, since we consider the LOSP mass to be 
broadly of order the weak scale, the LOSP freeze-out process cannot be ignored. Thus the nature 
of the signals depends on whether DM is dominantly produced via freeze-in of the FIMP or freeze- 
out of the LOSP and whether DM is the FIMP or the LOSP. This gives four scenarios which are 
summarised pictorially in Figure [3j 

Scenarios 1 and 2 (3 and 4) have FIMP (LOSP) dark matter, while cases 1 and 3 (2 and 4) have 
dark matter produced dominantly by freeze-in (freeze-out). Two schematic figures are shown for 
each of the four scenarios; the first gives the FIMP and LOSP spectrum, illustrating both freeze-in 
and freeze-out contributions to the dark matter. The figures on the right-hand side show the time 
evolution of the abundances of both the LOSP and the FIMP. Logarithmic scales are used so that 
the freeze-out and LOSP decay processes appear very sharp, while freeze-in occurs gradually. The 
masses of the particles involved in freeze-out and freeze-in have the same order of magnitude, so 
that freeze-out and the end of freeze-in occur at comparable times. At a later time, when the 
heavier of the LOSP and FIMP decays to the lighter, the abundance of the dark matter particle 
is boosted by the amount of the decaying particle abundance. 

The four scenarios are labelled according to the dominant DM production mechanism and the 
nature of DM. 
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Figure 3: Schematic representation of the four possible scenarios involving the freeze-in 
mechanism. The left-hand figures show the LOSP/FIMP spectrum with circles repre- 
senting cosmologically produced abundances. The large (small) circles represent the 
dominant (sub-dominant) mechanism for producing the dark matter relic abundance, 
a dotted (solid) circle signifies that the particle is unstable (stable), and a filled (open) 
circle corresponds to production by freeze- in (freeze-out). The right-hand figures show 
the LOSP and FIMP abundances as a function of time. The initial era has a thermal 
abundance of LOSPs and a growing FIMP abundance from freeze-in. The LOSP and 
FIMP are taken to have masses of the same order, so that FIMP freeze-in ends around 
the same time as LOSP freeze-out. Considerably later, the heavier of the LOSP and 
FIMP decays to the lighter. 

1. Freeze-in of FIMP DM The FIMP is the DM and the dominant contribution to the 
relic DM abundance is generated via the freeze-in mechanism. A small abundance of LOSP 
freezes-out which then decays late to FIMP dark matter. 

2. LOSP freeze-out and decay to FIMP DM The FIMP is again the DM but now the 

dominant contribution to the relic abundance is generated via the conventional freeze-out of 
the unstable LOSP which then decays to the FIMP. A sub-dominant component of FIMP 
DM arises from freeze-in. 

3. FIMP freeze-in and decay to LOSP DM The LOSP is the DM and the dominant 
contribution to the relic abundance comes from the freeze-in of a long lived FIMP which 
later decays to the LOSP. A sub-dominant component of DM arises from LOSP freeze-out. 

4. Freeze-out of LOSP DM The LOSP is again the DM but the dominant contribution to 
the relic abundance comes from conventional freeze-out of the LOSP. A small abundance of 
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FIMPs freezes-in and decays to give a sub-dominant component of LOSP DM. In the limit 
that this freeze-in abundance is small, the standard case of LOSP freeze-out is recovered. 

The accelerator and cosmological signals depend on the scenario. For example, production and 
decay of LOSPs at the LHC is possible in scenarios 1 and 2, while late decays during the MeV 
era of the early universe are possible in all scenarios. 



4.1 Long-lived LOSP decays at the LHC 

The hypothesis of dark matter generation via thermal freeze-out, e.g. WIMP dark matter, may 
be experimentally verified by a comparison of the observed dark matter density, i.e. Odm'* 2 , and 
properties of the WIMP measured at accelerators. The accelerator measurements may be used 
to infer the WIMP self-annihilation cross section owimp which, in the simplest cases, may be 
immediately employed to predict f2 W iMP^ 2 - 

Does a similar correlation exist for freeze-in? Freeze-in scenarios of dark matter typically 
contain both a stable dark matter particle and a very long-lived unstable state. In the case of 
FIMP dark matter, the unstable state is the LOSP. Since the LOSP can only decay dominantly 
via the same small coupling A that induces freeze-in, then there will be a relation between the 
LOSP properties, in particular its mass and lifetime, and the FIMP dark matter abundance. The 
crucial signature therefore arises from the accelerator production of the LOSP, followed by the 
observation of LOSP decay with a long lifetimd^ This signal is expected only for scenarios 1 and 
2; in scenarios 3 and 4 it is the FIMP that is unstable, but the FIMP production rate at the LHC 
will be negligible. 

The connection between the LOSP lifetime and the cosmological FIMP abundance is model 
dependent. For example, suppose that the LOSP is a scalar B and the FIMP freeze-in process 
occurs via a quartic scalar interaction leading to BB — > XX. In this case the LOSP decay 
must involve other small couplings that connect X to bath particles, so the relation between DM 
abundance and LOSP lifetime is lost. Here we concentrate on the particularly simple case that 
freeze-in occurs via the Yukawa coupling XBiB 2 X, where B\ is the LOSP and B 2 a lighter bath 
particle. The stability of X can be guaranteed by a Z 2 parity under which X and B\ are odd. 
The relic abundance of FIMP dark matter arises from LOSP decays, B\ — > B 2 X, and is computed 



in section 6.1 with the result of eq. (6.10). Imposing the WMAP constraint on the dark matter 



abundance, eq. (6.10) can be rearranged to yield a prediction for the LOSP lifetime 

„ 1ft _ 3 / m x \ /300GeV\V 10 2 \ 3/2 

t Bi = 7.7 x 10 3 sec<7 Bl — r 

VlOOGeW V m B 1 J \9*{m Bl )J 



(4.1) 



This result applies for scenario l, where the DM abundance is dominated by freeze-in. In scenario 
2 most DM arises from LOSP freeze-out, so that the fraction from freeze-in, /, is small, / C l. 
The collider signal from late decaying LOSPs persists, but now with a lifetime enhanced compared 



with eq. (4.1) by a factor l/f. 



investigations into the possible detection of long lived states at CMS are in progress, see for example 
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4.2 Reconstructed LOSP properties at the LHC 

As illustrated in Figure [3j regardless of the nature of DM (FIMP or LOSP) the LOSP will undergo 
conventional freeze-out. Consequently the dark matter abundance has two components, one from 
the freeze-in mechanism and the other from the freeze-out of the LOSP. The freeze-in contribution 
dominates only if the conventional thermal freeze-out cross section is large, av > 3 x 10 _26 cm 3 s -1 , 
corresponding to scenarios 1 and 3. 

Experiments at colliders may eventually measure the properties of the LOSP sufficiently well to 
determine that the annihilation cross section is indeed too large to give the observed dark matter 
from LOSP freeze-out, implying that the dominant contribution to the dark matter abundance is 
generated via an alternative mechanism such as freeze-in. 

It would be particularly interesting if LHC discovered prompt cascade decays of superpartners 
that terminated in a charged or coloured LOSP. Indeed, this would motivate a determined effort 
to search for a long lifetime associated with a LOSP, as in the previous sub-section. 

4.3 Enhanced direct and indirect detection signals of DM 

Consider the scenario where the LOSP is the dark matter and the freeze-in of a FIMP which 
later decays generates the correct relic abundance of LOSP dark matter (scenario 3 of Figure 
[3]). As discussed above this means that the conventional thermal freeze-out cross section for the 
LOSP DM is large: av > 3 x lCT 26 cm 3 s -1 . Hence, the indirect signals produced by LOSP DM 
annihilating via these large cross sections can be enhanced compared to the signals predicted for 
DM generated directly from freeze-out. 

Indirect detection experiments such as PAMELA [20], FERMI [21] and HESS [22] have recently 
reported deviations from background expectations of the proportion of positrons in cosmic rays. 
A common though perhaps premature interpretation of these experiments is the production of 
positrons resulting from WIMP annihilations. In order to match the observed fluxes one needs 
"boost factors" in the rates. Freeze-in provides a new avenue for building models containing 
boosted WIMP annihilation cross sections. 

We also note that increasing the interaction strength of the dark matter particles can lead to 
an increase in the scattering cross section relevant in direct detection experiments. This could 
provide further evidence in favour of the freeze-in mechanism if a signal was inferred in a region 
of parameter space not consistent with frozen-out dark matter. 

4.4 Cosmological decays during the MeV era &; perturbed BBN abun- 
dances 

As has been already discussed in some detail, freeze-in scenarios usually contain a metastable 
particle. This may be either the FIMP X itself or the LOSP Bi, whichever particle is the heavier 
one. It is well known that the existence of metastable particles in the early Universe may be 
constrained by the epoch of BBN. The light element synthesis of 2 H, 3 He, 4 He and 7 Li (as well as 
6 Li) may change drastically when hadronically and electro-magnetically interacting non-thermal 
particles are injected into the plasma [19], due to the decay of the metastable particles. Significant 
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deviations of the results of a standard BBN scenario occur only for decay times r > 0.3 sec. Com- 
paring this to eq. (4.1), and assuming a weak mass scale FIMP such that gf ~ g% ~ 100 as well 
as ~ 1, one arrives at a typical decay time of ~ 10~ 2 sec, implying no signatures from BBN. 
However, in case the FIMP particle couples to a large number of bath particles, i.e. gs 1 ~ 100, 
individual LOSP decay times may become ~ 1 sec, potentially perturbing BBN. 

Even larger deviations (i.e. increases by factor ~ 10 2 — 10 3 ) from the life time as given in 
eq. (4.1) may result if the assumed 2-body decays in eq. (4.1) are forbidden kinematically. As an 
example, consider the Yukawa interaction \XB\B 2 and a mass ordering > mj, such that 
the FIMP X is the dark matter. Assuming < ms 2 + mx any frozen-out LOSPs B\ may only 
decay via a 3-body decay since B 2 has to remain virtual. Three-body decay widths are suppressed 
compared to 2-body decay widths resulting in longer decay times. Given this, the typical decay 
time could move into the range tb 1 ^ 3 — 3000 sec, the regime where BBN may be significantly 
perturbed. 



4.5 Generating a warm DM component: Erasure of small scale density 
perturbations 

Decay produced particle DM is often warm/hot, i.e. is endowed with primordial free-streaming 
velocities leading to the early erasure of perturbations, due to the kinetic energy imparted on the 
decay product during the decay itself. As at least part of the dark matter in freeze-in scenarios is 
produced by the late decay of the LOSP to the FIMP, or vice versa, and since neither LOSPs nor 
FIMPs may thermalize below cosmic temperatures T< 1 — lOMeV, freeze-in scenarios may come 
in the flavour of warm or mixed (i.e cold and warm) DM scenarios. How warm depends strongly 
on the decay time, and the mass ratio of mother and daughter particle. 



5 Abundance "Phase Diagrams" 

In this paper we argue that freeze-in can provide a possible alternative thermal production mech- 
anism to the much studied freeze-out mechanism. We now sketch "abundance phase diagrams" 
that show regions of mass and coupling where each mechanism dominates the production of the 
relic abundance. We find that the following simple framework allows four distinct physical be- 
haviours to account for the dominant production of dark matter. There is a thermal bath that 
contains species that are kept thermally coupled during the era of interest by interactions with 
strength of order unity. In addition there is a particle, A, that couples to the bath particles by a 
renormalisable interaction with a coupling A. 

We assume that there is an unbroken symmetry that leads to the stability of the dark matter 
particle, which may be A or a bath particle. This symmetry is carried by some of the bath 
particles, and possibly also by X. As the coupling between X and the bath decreases, late decays 
often become relevant for the production of dark matter: decay of any X that is produced when 
X is not the dark matter, or decays to A if A is the dark matter. 

The interaction between A and the bath has three possible forms: a quartic or trilinear scalar 
interaction, or a Yukawa interaction. Furthermore the final dark matter abundance depends on 
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Figure 4: In (a) we show the contours of Qh 2 as a function of the mass mx and 
coupling A for the case of a quartic interaction. The plane can be divided into two 
"phases": X freeze-out, phase (I), occurs for large coupling and X freeze-in, phase 
(III), occurs for weak coupling. In (b) we take a slice at fixed m x ~ v and plot the 
variation of flh 2 as a function of the coupling A. 



whether X is a fermion or boson, and on how many X fields appear at the interaction. For each 
such scheme one can compute the dark matter abundance in terms of the coupling and the masses 
of the relevant particles, and we find that there are just four possible behaviours 

• (I) Freeze-out of X. 

• (II) X decouples with a full "T 3 " number density. 

• (III) Freeze-in of X. 

• (IV) Freeze-out of bath particle. 

Hence, the parameter space of any such theory can be split into four "phases" according to which 
behaviour yields the dominant dark matter abundance. Two of these phases correspond to the 
well-known case of freeze-out (I and IV). However, the introduction of very small renormalisable 
couplings to a theory introduces the possibility of two new phases; one with freeze-in (III) and 
the other having decoupling with a full relativistic abundance (II). 

The above list is ordered so that, as the interaction strength between X and the bath is 
decreased, one passes down the list, although all four "phases" do not appear in all theories. 
Below we give two illustrative examples of this "phase" behaviour, then we argue that for very 
small values of A the "phase" structure takes a universal form. 

5.1 Phase diagram for a quartic interaction 

Suppose that X is a scalar and the interaction with the bath is via a quartic coupling to a bath 
scalar B 

L Q = \B ] BX ] X. (5.1) 
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The stabilising symmetry is a Z 2 with X odd and B even. X is the lightest particle with odd 
charge, and hence is stable and yields the dark matteij^J B is unstable, decaying rapidly to 
other bath particles. X has no Standard Model gauge interactions, and has a cosmological relic 



abundance determined by the interaction eq. (5.1) 



For simplicity we consider only the case rax > m B, so that the relic X abundance depends 
on (mjjf,A). In Figure [4^, contours of Vth 2 are shown in this plane. The plane can be divided 
into two "phases": X freeze-out, phase (I) occurs for large coupling, A 2 > m x /M Ph above and 
to the left of the diagonal line; while X freeze-in, phase (III), occurs for weak coupling with 
A 2 < mx/Mpi. These are the only two behaviours, there are no regions with phases (II) and (IV), 
and the plane is equally divided in logarithmic space between them. The parametric equations 



for the abundances in these regions are given by eqns. (2.4) and (2.2), with m = mj. The phase 
boundary occurs when Y FO ~ Y FI ~ 1: for freeze-out, A is sufficiently small that freeze-out occurs 
right at Tpo ~ mx and not below, while for freeze-in A is sufficiently large that a full thermal 
abundance just freezes-in by Tpi ~ mx- In both phases, Y drops with distance from the boundary. 
In the freeze-out region contours of fixed Qh 2 are straight lines: A 2 oc (1/Qh 2 )(mx/Mpi) 2 , while in 
the freeze-in region these contours are independent of mx'- A 2 oc Qh 2 . Throughout the entire plane 
the final abundance is set in the radiation dominated era; although along the phase boundary the 
abundance is being set during the transition from radiation to matter domination. 

The observed value of flh 2 is shown by a contour with a dotted line. Along this contour, the 
transition from freeze-out to freeze-in occurs for a particle of mass v 2 /M u just freezing-in a thermal 
abundance and immediately freezing-out at a temperature v 2 /M u , where M u ~ 10 16 GeV is the 
scale of gauge coupling unification. This special case gives hot dark matter and is excluded; it is 
hard to engineer as it would require X to be interacting with electrons, photons or neutrinos at 
the eV era. Moving away from this special point on the contour, the dark matter mass increases, 
becoming warm and then cold dark matter, whether on the freeze-out or freeze-in side. 

There are two particularly interesting parts of this contour: the first has A ~ 1 and mx ~ v, 
corresponding to the well-known WIMP case where the physics of dark matter depends on a single 
scale, v. The second has A ~ v/M u and is independent of rax- This is the freeze-in region, and 
within it, the case of mx ~ v is particularly interesting since in this case there are only two 
scales associated with the dark matter, the weak and unified scales. These two special cases, both 
having mx ~ v, are shown by solid dots. One has A ~ 1 and corresponds to the WIMP case, 
while the other has A ~ 10 -12 ~ v/M u , and is the FIMP case. A slice through the parameter 
space with mx ~ v is shown in Figure [4]d, which plots the prediction for Vth 2 as a function of A 2 . 
While the WIMP case is presumably related to the physics of electroweak symmetry breaking, 
FIMP freeze- in with rax ~ v and A ~ v/M u may be related to the moduli of extra dimensional 
theories at the unified scale, such as string theory. Providing the initial abundance of the moduli 
is suppressed, for example by well-known solutions to the moduli problem, a modulus or modulino 
can be frozen-in during the electroweak era to form dark matter. For both WIMP freeze-out and 
moduli freeze-in, the amount of dark matter is parametrically given by 



3 The phenomenology of this interaction has been previously considered in ref. |23j . 
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Figure 5: In (a) we show the abundance phase diagram with contours of Qh 2 as a 
func tion of the coupling A and the mass mx for the case of a Yukawa interaction, 



eq. (5.3). If A 2 > yJmx/Mpi, X undergoes conventional freeze-out giving region I; 



while in region II, with mx/mpi < A 2 < \prnxjMp~h X decouples from the bath giving 
a yield Yx ~ 1. In region III, with mx/mpi > A 2 > (mx/Mpi) 2 , Yx never reaches unity 
and freeze-in provides the dominant contribution to dark matter. For A 2 < (mx/Mpi) 2 
the dominant mechanism generating dark matter arises from ipi freezing out and then 
decaying to X + ip 2 , giving region IV. In (b) we take a slice at mx = 1 TeV and plot 
the variation of Vth~ function of the coupling A. 



where no distinction is made here between M u and Mp\. 



5.2 Phase diagram for a Yukawa interaction 

Suppose that the scalar, X, interacts with the bath via a Yukawa interaction with two fermions 

L Y = XfaifaX, (5.3) 

where m-i <C mx < m\, with mx and mi of the same order of magnitude. The stabilising 
symmetry is a Z 2 carried by both X and ipi, which is the lightest bath particle odd under the 
symmetry. Hence, the dark matter is X and it has an abundance with an order of magnitude 
determined by (mj,A). Contours of Vlh 2 are shown in Figure^, and all four phases discussed 
above occur, in the regions labelled (I), (II), (III) and (IV). 

The initial X abundance is assumed to be negligible, and as we progress from phase I to 
phase IV the freeze-in process becomes successively less efficient. If A 2 > mx/Mpi, corresponding 
to regions I and II, a full thermal X abundance is produced at temperatures above mj. If 
A 2 > \fm~xjMph X undergoes conventional freeze-out via the reaction XX — > tp2^>2, giving 
region I; while in region II, with m x /M P i < A 2 < ^Jm x /M Ph X decouples from the bath at 
T > mx, giving a yield Yx ~ 1. On the other hand, with A 2 < mx/Mpi, the freeze-in process, 
'ipi ~ * ' l l J 2X, is less efficient, so that Yx never reaches unity in regions III and IV. In region III, 
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Figure 6: Fixing the bath particle mass m x = 1 TeV the phase diagram shows more 
structure. For mx < mi, X is the dark matter as in Figure |5j A new feature compared 
to Figure [5] emerges for mx > mi, where ipi is the dark matter. In this case dark matter 
can be dominantly produced by the freeze-in of X at temperatures T > 1 TeV, and the 
subsequent decay of X — > ipi + ip 2 . Here the decay has to occur after the freeze-out of 
■ipi, otherwise the ipi contribution due to X-decay will be reprocessed. The X-decay 
contribution of ipi dark matter dominates the freeze-out contribution in the triangle 
defined by the boundaries of zone III and the condition mx > m\. A slice through this 
figure with mx fixed at a TeV gives the same dependence as shown in Figure [5}d. 



with mx/Mpi > A 2 > (mx/Mpi) 2 , freeze-in provides the dominant contribution to dark matter, 
while, for A 2 < (mx/Mpi) 2 the dominant contribution to dark matter arises from ipi freeze-out, 
ipiipi — > lighter bath particles, followed by ipi — > Xip2, giving region IV. 

Figure [5Jd shows a slice through the two dimensional parameter space of Figure [5^, with mx 
fixed to 100 GeV. For a large range of A too much X dark matter is produced. But, as with the 
quartic interaction, there are two very interesting values of the coupling which yield the observed 
abundance of dark matter, WIMPs with A ~ 1 and FIMPs with A ~ 10~ 12 . For these two special 
cases VLh? passes through the observed value, as shown in by solid dots Figure [5^,. FIMPs with 
a weak-scale mass lie close to the phase boundary between X freeze- in and ipi freeze-out. This 
important case, motivated, for example, by the possibility that X is a modulus, will be examined 
in some detail in the next two sections. Many of the observational consequences of weak-scale 
FIMPs arise because there are two broadly comparable production mechanisms, involving both 
freeze-in and late decays. 

In Figure [5j we have mx and m\ of the same order of magnitude as we scan over different 
values for mx- We can develop this scenario by fixing mi = 1 TeV while still allowing mx to 
vary. The resulting phase diagram is shown in Figure [6j Once again the diagram is split into four 
different regions but now these regions, in particular region IV, have more structure. The first 
feature to notice is the vertical dividing line at mx = mi = 1 TeV which separates the two cases 
where mx is the lightest and therefore the dark matter (left hand side) and mi is the lightest and 
therefore the dark matter (right hand side). 
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Figure 7: Independent of the interaction, there is universal behaviour below the upper 
shadowed line, corresponding to small coupling, A 2 <C mx/Mpi. For such couplings 
the dominant production of dark matter arises from freeze-in of X if A > mx/Mpi, 
and freeze-out of a bath particle if A < m x /M P i. 

As before, region I corresponds to the conventional freeze-out of X via the reaction XX — > ^2^2 
and region II corresponds to where the X particle decouples from the thermal bath with a yield 
Y x ~ 1. Both these regions are cut off at m x = mi due to the fact that when m x > mi, X is no 
longer stable and decays to bath particles. Consequently, the dark matter abundance arises from 
ipi freeze-out, V'iV'i — > lighter bath particles. 

Moving down the diagram to smaller values of the coupling A, region III is reached and the 
freeze-in mechanism dominates. For m x < m 1; X is the dark matter with the freeze-in process, 
■01 ~~ * ^2^-1 generating the dominant contribution to the relic abundance. Moving into the region 
where m x > mi we now freeze-in an abundance of X particles via the process ipi +1P2 —> X. The 
resulting abundance of X particles then decays back to ipi which forms the dark matter. Within 
this region as we move to larger values of m x at constant A we decrease the abundance of X 
particles frozen-in but we also decrease the lifetime of the X particles. Eventually, the lifetime 
becomes so short that the X particles will decay back to the bath particles before ipi freezes out. 
This means that the relic abundance of ipi particles is once again determined by the freeze-out of 

Moving to even smaller A we move into region IV where the freeze-in mechanism becomes less 
efficient. For m x < mi the relic abundance of X particles is determined by the freeze-out density 
of tpi particles which then decay to X particles after freezing out. For m x > m 1 the X particles 
play no role in determining the relic abundance of dark matter. The ipi particle is the dark matter 
with it's relic abundance determined by conventional freeze-out. 
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5.3 A universal phase diagram at small coupling 



In Figure [5^, the two phases at lowest coupling correspond to X freeze-in, III, and ipi freeze-out, 
IV. We find that this phase structure at small coupling results in a very wide class of theories. 
Consider theories where X is either a scalar, A x , or a fermion, ip x , with linear couplings to bath 
scalars (A) and fermions (ip) that are renormalisabk^J so that the possible interactions are 

XA x tpijj, \m x A x A 2 , XA X A 3 , and \ipxi>A. (5.4) 

The stabilizing symmetry is Z 2 , and A x ,i/j x may be even or odd (although, if it is even, it 
must have a significant fraction of its coupling to bath states that are odd). Any choice of parity 



may be made for A and ip as long as the interactions of eq. (5.4) are symmetric. We assume 



that Z2-odd bath particles in the interactions of eq. (5.4), as well as the lightest bath particle 



odd under the Z 2 , B, have masses of the same order of magnitude as m x , and that Z 2 -even bath 



particles in the interactions of eq. (5.4) do not have masses parametrically larger than m x . 

For this wide class of theories, the phase diagram is shown in Figure |7[ for A 2 < mx/Mpi. 
The dominant production of dark matter, whether X or bath, arises from freeze-in of X, if 
A > m x /Mpi, and freeze-out of a bath particle, if A < m x jMp\. This generality follows from 
two results. Since B, the lightest Z 2 -odd bath particle, has order unity couplings to lighter bath 
particles, its freeze-out yield is always Yp ~ itib/Mpi ~ m x /Mpi. Freeze-in of X via 2^2, 



1 — > 2 or 2 — > 1, involving any quartic, trilinear or Yukawa interaction of eq. (5.4), leads to a yield 
Y x ~ X 2 M Pl /m x . 

^From Figure [7] we see that this class of theories contains the interesting weak-scale FIMP 
case, m x ~ v and A ~ v/M u . This leads to a generically realistic abundance of dark matter, 
arising from two broadly comparable mechanisms — X freeze-in and B freeze-out — as shown by 
the FIMP label on the phase boundary in Figure [7| 



6 Freeze-In Calculation 

We now turn to the calculation of the frozen-in dark matter density. The freeze-in process is 
dominated by decays or inverse decays of bath particles to the FIMP depending on whether or 
not the FIMP is the lightest particle carrying the conserved quantum number that stabilises the 
DM. Processes such as 2 — > 2 and other scatterings give sub-dominant contributions to the frozen- 
in abundance in almost all cases of interest as the diagrams are suppressed by both additional SM 
gauge or Yukawa couplings and numerical factors of order few x 7r 3 arising from the additional 
phase space integrals. 

6.1 Direct freeze- in of dark matter 

Consider the case where the FIMP X is the dark matter particle itself. If X has a coupling to two 
bath particles, XXBiB 2 , then, for mp 1 > mp 2 +m x , the dominant freeze-in process is via decays 

4 For simplicity we ignore the quadratic interactions AAxipx*^, ^A X A 2 , XAipxtpx, and \mxA\A. A similar 
analysis holds in the more general case. 
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of the heavier bath particle: 

B x -> B 2 X. (6.1) 

At temperatures much higher than m^, the yield during a Hubble doubling time at the era of 
temperature T scales as 

y^(t) cc MnBsI*, (6 . 2) 

so that it is strongly IR dominated, and shuts off only once the temperature drops below M Bl . 

We can be more precise by solving the Boltzmann equation for the number density of X 
particles in this case: 



h x + 3Hn x = J dUxdU Bl dU B2 (2n) 4 5 4 (px + Pb 2 ~ PbJ 

x [\M\l^ B2+x f Bl {l±f B2 ){l±fx)-\M\l 2+x ^ B J B J x {l±f Bl )], (6.3) 

where dUi = d 3 pi/(27r) 3 2Ei are phase space elements, and fi is phase space density of particle i, 
such that 

* = A / #Pfi (6.4) 



2vr 3 „ 

is the total particle density of species i possessing g»j internal spin degrees of freedom. It is 



implicitly assumed that in eq. (6.3), as well as eqns. (6.15) and (6.16) below, the squares of 
the matrix elements are summed over final and initial spin of the participating particles without 
averaging over the initial spin degrees of freedom. We assume the initial X abundance is negligible 



so that we may set fx = 0, such that we may neglect the second term in eq. (6.3). Using the 
definition for the partial decay width T Bl of B\ — > B% X and neglecting Pauli-blocking/ stimulated 
emission effects, i.e. approximating (1 ± f Bi ) pa 1, we can rewrite the Boltzmann equation as 

h x + 3Hn x « 2g Bl [ d!i Bl T Bl m B J Bl = g Bl [ t/fo ^ Bl ^ Bl (6.5) 
J J (2viT 7 Bl 

where j Bl = E Bl /m Bl . The bath particles are assumed to be in thermal equilibrium and so 
approximating f Bl = (e EB i/ T ± 1) _1 by e~ EB \l T and converting the integral over momentum space 
into an integral over energy we have 

I q TT ^ f ^P Bl f Bl T Bl f°° m Bl T B 2 2 1/2 -EgjTjrp 

n x + Zn x H w g Bl J ^ — = g Bl —^-(E Bl - m B J 1 e ^ dE Bl 

= gBl " 2 f Bl TK 1 (rn Bl /T). (6.6) 

where K\ is the first modified Bessel Function of the 2nd kind. Rewriting in terms of the yield, 
Y = n/S and using T pa —HT, applicable when the variation of total plasma statistical degrees 
of freedom with temperature dg/dT pa approximately vanishes, we have 



x 



Tmax 9 Bl mlT Bl K^msjT) 

SH— dT > (6 ' 7) 
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with S = 27r 2 gfT 3 /45 and H = 1.66yg£T 2 /Mpj. Using x = m/T we can rewrite the integral as 

lax 

K 1 {x)x^dx 1 (6.8) 



Y „ 45 g Bl M Pl T Bl 
X ~ (1.66)4tt* m^VS 

where Mpi is the (non-reduced) Planck mass. Doing the x integral with x max — oo and x m i n — 
we finally arrive at the result 

Y^ 2 ~ 135 ^ = f , (6.9) 

8tt3(1.66)^v^ V rn 2 Bi J 

where ,p are the effective numbers of degrees of freedom in the bath at the freeze-in temperature 
T ~ m Bl for the entropy, S, and energy density, p, respectively. Thus the abundance depends on 
the three quantities mx,mB 1 and Tb 1 and has the form 

^^1.09X10^^, (610) 

gWa* m B 1 

This is the density of X produced by a single bath particle species. In full FIMP models it is quite 
likely that a number of bath particles could have similar interactions with the FIMP particle. 
For example, if X arises as an i? p -odd modulino state from the SUSY-breaking sector of a string 
theory, then as discussed in Section [3j X is likely to couple to most if not all of the MSSM 
spectrum. In this case further contributions to the relic abundance are generated with the same 
dependence on the bath particle's partial decay width and mass, giving 



n x h 2 \ to t « - ~j = j> m x 2^ ( 6 - n ' 



and the overall interaction strength will need to be smaller in order to generate the correct abun- 
dance of FIMP dark matter. If we approximate this enhancement by an effective number of active 
degrees of freedom, gbath, each with mass mg, then the required interaction strength becomes 



A.1.5xl0-feV / 7^Y ,/4 (^)- 1/2 



mx) \ 10 2 J V 10 2 
6.2 Decays of frozen-in FIMPs to dark matter 



(6.12) 



The alternative possibility is that the particle that freezes-in is unstable and decays to dark matter. 
Here we study the simplest possibility that the interaction responsible for freeze-in also yields the 
decay. If X has a coupling to two bath particles, XXB 1 B 2 , then, for m x > m Bl + mB 2 , the 
dominant freeze-in process is via inverse decays of X: 

B X B 2 -> X. (6.13) 

At temperatures much higher than mx, the yield during a Hubble doubling time at the era of 
temperature T scales as 

v (rrt\ M Pi m x^x fR , A , 
Y 2 _ >1 {T) oc — , (6.14) 
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and is again strongly IR dominated, shutting off as the temperature falls below m x . 

Assuming once again that the initial abundance of X particles is zero and therefore setting 
fx — the Boltzmann equation for this process can be written as 



h x + 3Hn x ~ J dIl x dIl Bl dIl B2 (2iT) 4 5 A (px ~Vb y ~ Pb 2 ) \ m \b 1+ b 2 ^x /bi/b 2 - 



(6.15) 



Assuming CP invariance we may set \M\ 2 Bi+B2 ^ x = \M\ x _ tBi+B2 and by invoking the principle of 
detailed balance we can rewrite the Boltzmann equation as 

h x + Sn x H « J dU x dU Bl dU B2 (27r) 4 5\px - p Bl - Pb 2 ) \M\ x _ Bl+Ba f x q , (6.16) 

where f x is the X equilibrium phase space distribution approximated again by f x ~ e~ Ex ^ T . 
Comparing eqns. (6.16) with (6.3) we can immediately write down the resulting form of the X 
yield as 

135 fMpiT x ' 



•2-*l 



(6.17) 



8tt 3 (1.66)^v^ V m x 

where T x is the partial width of X — ► B\B2- 

Assuming that B\ is the DM particle, and in addition that the freeze-in contribution coming 
from decays of X dominates the conventional freeze-out abundance of Bi, the final DM density is 

a Bl h>» im *}?' (6.i8) 

Here, crucially, we have assumed the decay of X — > BiB 2 occurs at a time after the freeze-out 
of B\ so that the density does not get reprocessed. Taking T x = X 2 m x /8n, the required dark 
matter density occurs for a coupling of size 



Although eq. (6.18) is very similar in form to eq. (6.10), as was to be expected, the physics is 
quite different. 



6.3 Freeze-in by 2^2 scattering 

Finally we present the calculation of the FIMP relic abundance in the case where the FIMP, X, 
is a scalar and interacts with three scalar bath particles Bi, B 2 and B 3 via the operator 

^4-scalar = XXBiB 2 B 3 . (6.20) 

Considering this interaction we can calculate the resulting FIMP yield using the Boltzmann equa- 
tion 

h x + 3n x H w 3 / dU Bl dU B2 dU B3 dU x (27r) 4 5 4 (p Bl +p B2 -p B . i -p x )\M\ BiB2 ^ B3X f B J B2 , (6.21) 
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where the factor of 3 accounts for the fact that we can have B\B 2 — > B 3 X, BiB 3 — > i?2-^ and 
— ► contributing to the FIMP yield all with the same rate. We assume that the masses of 
Bi, B2 and B3 are negligible compared to the FIMP particle mass. We can rewrite the Boltzmann 
equation as a one dimensional integral, [21] , 

noi roc 

n x + 3n x H^—: / ^ ds dQP BlB2 P BaX \M\ 2 BiB ^ BsX K^/T)/^, (6.22) 

Jm? x 

where s is the centre of mass energy of the interaction at a temperature T and 

__ [s - (rm + m,) 2 ] 1/2 [ g ~ K ~ ^)T /2 (a ~\ 



The matrix element is \M\ B B _^ B x = X 2 leaving 

3TA 2 f°° - - 

h x + 3n x H^—— dsis-m^K^/T)/^, (6.24) 

Doing this s integral and using the definition for the yield, Y = n/S we have 
dY x 3X 2 T 2 m x K^mx/T) _ 3X 2 K 1 (m x /T) 45M pl m x 



dT SH 128tt5 L66 T*gSy/tf 256tt 7 



(6.25) 



Changing variables from T to x = m x /T and again doing the integral x integral under the 
approximations that x max = 00 and x m i n = we finally arrive at the result 



■x 



l35X2M « r xKi(x)dx= . (6 . 26) 



256n 7 gSy/gP(1.66)m x Jo 512ttVV gS(1.66)m x 

The relic density of X FIMPs is then given by 



, 2m x Y x 1.01 x 10 24 _o 



Finally to generate the required relic abundance we need 

3/4 

10 2 , 

larger than the corresponding value for the three body interactions. 



A~lx 10- 11 ( ) , (6.28) 



7 Comments and Discussion 

Having outlined the elementary theory of the freeze-in mechanism and FIMP phenomenology we 
now address some of the complications and issues that can arise in complete theories or ones where 
the FIMP sector is not just one particle. 



24 



7.1 Higher dimension operators and FIMPs from GUTs 

In addition to the example candidates presented in section [3] an important application of the 
freeze-in mechanism is to non-renormalisable higher dimensional operators (HDO)s containing at 
least one Higgs-like state. Consider the operator 

C-HDO = —(fl^-fnjX^!^, (7.1) 

where a is an 0(1) coupling, X is our FIMP state, ipi and ip2 are two fermionic bath states and 
the ipiS are Higgs-like states, also assumed to be in thermal equilibrium, that will gain vacuum 
expectation values (VEVs) at some energy scale below the high scale M. These higher dimension 
1/M suppressed operators involving Higgs VEVs (e.g. for n = 1 electroweak scale VEVs or for 
n > 1 larger-scale VEV^]) can lead to interactions similar to the dimension-4 Yukawa interaction 



of eq. (5.3) with the appropriate coupling size to generate the correct freeze-in abundance of dark 
matter. 

For example, dimension five 1/Mcc/T-suppressed operators involving the SM Higgs VEV can 
lead to interactions with SM-singlets contained in GUT representations giving FIMPs with inter- 
action strength, A ~ v/M GUT ~ 10~ 13 . Alternatively, dimension six l/M^^-suppressed operators 
involving two intermediate scale VEVs can also give rise to suitably sized FIMP coupling^} How- 
ever if the excitations of the Higgs-like states are not super-massive and are in thermal equilibrium 
there is a UV contribution to the dark matter yield that limits the applicability of the IR freeze-in 
mechanism. 

For the sake of clarity we specialise to the case of n = 1 with one Higgs state gaining an 
electroweak scale VEV and mass. Expanding cp\ around this VEV we have the following operators 

(XV Ct 

£hdo = jjXfafo + —¥1X^2- (7.2) 

Under the assumption that > m^ 2 + mx the IR dominated freeze-in yield is determined by 
the decay ip\ — *> 1P2X with rate determined by the coupling a' = av/M. The resulting yield can 



be directly read from eq. (6.9) and has the form 



_> 



where we have approximated the rate as T t p 1 « a' 2 This contribution is IR dominated 
by temperatures close to the mass m^. However, we also get a contribution to the yield of X 



coming from the non-renormalisable interaction in eq. (7.2) via a four particle interaction such 
as ipi^pi — > Xip2 and the important point is that this contribution is UV dominated, so this 
contribution will depend on unknown UV physics such as the reheat temperature Tr. 



5 For n = 0, the resulting operator is renormalisable and has an 0(1) coupling. This interaction will not lead to 
freeze-in due to the large coupling and so we assume it is absent due to symmetries. 

6 For existing models employing 1 /Mqut operators in the context of decaying dark matter see for example [25 . 
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What constraint is there on T R for the yield to be determined primarily by IR physics? Consid- 
ering the HDO interaction only we can calculate the resulting yield using the Boltzmann equation 

n x + 3n x H « J dU^dU^dU^dUx^TiYS 2 ^ + p v - p^ 2 - p x ) \M\\ i(p ^ 2X f v f^. (7.4) 

Manipulating this equation we can write 

n x + 3n x H* j dsdn^M^^K^s/T), (7.5) 

where s is the centre of mass energy of the interaction at a temperature T and we have approxi- 
mated the masses of the relevant particles to be negligible compared to the temperature at which 
we are working. In this limit the matrix element is \M\*,^ 2X = a 2 /M 2 s leaving 

Try 2 r°° 

hx + 3nxH ~ 512^1 dss3/2K ^/T), (7.6) 



Doing this final integral and using the definition for the yield, Y = n/S we have 

dY uv 1 TV 



dT SHT 16tt 5 M 2 ' 

Now performing the final T integral which is dominated at the highest temperature, T R . 



(7.7) 



^ 0AT R a 2 M, 

Thus in order for the IR contribution to dominate over that arising from the UV we need 

TT~ ^ TfT > 1, (7-9) 

translating into an upper bound on the reheat temperature given by 

^tt 3 7; 2 

T R <— . (7.10) 

With = 150GeV, the maximum value for T R is around 20 TeV giving a potentially serious 
restriction on this type of model. 

The examples we have considered in the main body of this paper avoid this problem as the 
interactions determining the freeze-in dynamics come from renormalisable or super-renormalisbale 
operators. Moreover, if the excitations of the "Higgs" field that converts the HDO into a d < 4 
operator approach the intermediate scale, as might be expected if the associated VEVs are at 
this scale, then no significant restriction on T R arises. Overall, there are many examples of non- 
renormalisable operators and their significance in each case is somewhat model dependent. 

It needs to be emphasised that similar issues occur with freeze-out in the presence of non- 
renormalisable operators involving a thermally decoupled particle that is either stable or long 
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enough lived that it decays after the freeze-out of the would be dark matter. Since extensions 
of the SM involving GUT-unification, gravity, or extra dimensions often have such particles and 
couplings, freeze-out theories of dark matter genesis similarly suffer from restrictions on Tr when 
viewed in this broader context. The most well studied example involves the gravitino where 
the yield Yq ~ 10~ 12 (TR/10 10 GeV) also approximately increases linearly with Tr, leading to 
important limits on the reheat temperature for weak-scale gravitino mass, though there are also 
other examples of particles interacting via non-renormalisable operators with analogous limits on 
Tr, see e.g. [26]. 

7.2 Freezing-in a FIMP sector 

Throughout this paper we have focussed on the case of one FIMP species. We note here some of 
the potential extensions of this idea by considering more complicated FIMP sectors. For example, 
consider the case where we have one FIMP species, X, coupled to a thermal bath via some 
renormalisable operator. It is through this operator that an abundance of X is frozen in. In 
addition to this suppose X has further interactions with some other states, X it which have even 
weaker or no direct interactions with the thermal bath. The strength of the interactions between 
the XjS and X could be moderate or even large. Consequently, once an abundance of Xs are frozen- 
in these extra interactions could bring the X^s and X particles into partial thermal equilibrium 
and within this sector a secondary freeze-out process could take place leaving a relic abundance of 
the lightest state in this sector. If the lightest state is stable it can form the dark matter but if it is 
not it decays back into an LSP or similar in the visible sector. This process of secondary freeze-out 
can substantially alter the final dark matter abundance, changing the observable phenomenology 
associated to the FIMP. 

A further variation arises when there is more than one FIMP particle, for example we can 
imagine that there are several moduli (perhaps with different masses) that couple via some operator 
to particles in the thermal bath. The couplings of these moduli are feeble and so each could have an 
abundance frozen-in. If the moduli have different coupling strengths (or masses) the abundances 
of each moduli will be different. Depending on the mass spectrum these frozen-in moduli could 
be the dark matter or decay to the real dark matter particle. 

8 Conclusions 

The nature and origin of the dark matter in the universe is unknown. While there are many 
theories of DM, there are rather few cosmological production mechanisms, and even fewer that 
can be subjected to precision tests. For example, the non-thermal coherent oscillation of a scalar 
field, such as an axion, always involves an unknown initial field amplitude. Thermal mechanisms 
that are independent of initial conditions are particularly interesting, and the decoupling of a 
heavy particle species from the thermal bath has received enormous attention. Such decoupling 
can occur when the particle is relativistic, as with the neutrinos, or non-relativisitic, as with 
baryons. The former leads to hot dark matter, so that the latter case of thermal freeze-out has 
been widely studied for DM. 
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At first sight these appear to be the only two ways of thermally producing DM without any 
sensitivity to initial conditions. If a massive particle of the thermal bath is to survive the early 
universe with a significant abundance it must decouple, and this can either happen when the 
particle is relativistic or non-relativistic. We have argued that there is one other possibility: the 
massive particle may have a negligible initial abundance and may be produced by collisions or 
decays of particles in the thermal bath. This is only a new possibility if the massive particle 
does not reach thermal equilibrium, and hence requires that it is a FIMP, interacting very feebly 
with the bath. Since the FIMP abundance is heading towards thermal equilibrium we call this 
production mechanism "freeze- in." 

Freeze-in requires a very small coupling A between the FIMP and the bath 



A 



2 T eq itifi 



Mpi m DM 



where T eq is the temperature of matter-radiation equality, mpi is the mass of the heaviest particle 
involved in the freeze- in reaction and ttldm is the mass of the DM particle. The mechanism 
works for a very wide range of FIMP masses: to be produced it must be lighter than the reheat 
temperature after inflation, and for DM to be cold it must be heavier than a keV. Indeed the 
range of theories involving FIMP freeze-in is so large that in this paper we have concentrated on 
FIMPs with a mass of order the weak scale, such as moduli or right-handed sneutrinos. 

It is not surprising that freeze-out is so popular: it presumably occurred with baryons (al- 
though with the added complication of a chemical potential) and it works with dimensionless 
coupling parameters of order unity. Indeed, with couplings of order unity, if the only mass scale 
in the annihilation cross section is that of the DM particle, then this mass is predicted to be of 
order the weak scale. On the other hand freeze-in of a FIMP requires a special situation. The 



freeze-in reaction must involve a small dimensionless coupling parameter, eq. (8.1); furthermore, 
UV sensitivity may reappear unless higher dimensional operators are suppressed or the reheat 
temperature after inflation is quite low. However, small couplings for both renormalisable and 
higher dimensional operators occur very easily, for example by approximate symmetries or small 
wavefunction overlaps in higher dimensions. Most importantly, precisely because FIMPs neces- 



sarily have very small couplings to the bath, eq (8.1), they offer the prospect of exotic signals of 
dark matter generation. 

If the FIMP is the lightest particle carrying the DM stabilising symmetry it is the DM, so that 
the lightest observable sector particle (LOSP) with this symmetry is unstable and decays to the 
FIMP with a rate suppressed by A 2 . Alternatively it may be the LOSP that is DM, and in this case 
it is the FIMP that has a very long lifetime. The signals associated with the freeze-in mechanism 
all revolve around the long lifetime of decays between the FIMP and LOSP. In particular, we 
explored the signals that arise when the FIMP, LOSP and freeze-in masses are all of order the 
weak scale. 

In many theories a small coupling may arise from a parametric form that is linear in the weak 
scale v , A ~ v/Mpi. In fact, in many models the Planck scale will be replaced by the string or 
compactification scale, as occurs for certain moduli, providing a better explanation for the order of 
magnitude of A. The temperature of matter-radiation equality from DM produced by the freeze- in 
mechanism then has the same parametric form as in the WIMP case: T eq ~ v 2 ' jMp\. This not only 
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explains why freeze-in can yield the observed abundance, but suggests that two mechanisms for 
DM production - FIMP freeze-in and LOSP freeze-out - yield broadly comparable abundances. 

The cosmological signals of the FIMP arise from the late decays between the FIMP and the 
LOSP and, with masses of order the weak scale, the parametric form of the lifetime in the simplest 
theories is 

Mpi . , 

r ~ — — ~ 1 sec. (8.2) 
T eg v 1 ; 

A precise computation shows that if the FIMP is frozen-in via two body decays or inverse decays 



that involve the LOSP, then the lifetime is closer to 10~ 2 seconds, as shown in eq. (4.1 ). However, 
this should be viewed as a lower bound on the lifetime of decays between the FIMP and LOSP. 
Further suppression of the decay rate occurs if the freeze-in process involves bath particles other 
than the LOSP or if the decays are to three or more particles. Thus the lifetime may well be in the 
region where the decays can alter the nuclear abundances produced during BBN, possibly solving 
the 7 Li problem, and may even be long enough to allow a component of the DM to be warm. 

If FIMP freeze-in dominates over LOSP freeze-out, then the LOSP annihilation cross section 
is larger than with conventional freeze-out. This leads to an important astrophysical signal of 
LOSP DM: the DM annihilation in the galactic halo is boosted relative to the conventional case, 
leading to enhanced rates for indirect detection signals of photons, leptons and anti-protons. 

Finally, FIMP dark matter yields important new collider signals. The exotic particles of the 
LOSP sector can be produced at the LHC and will rapidly decay to the LOSP, which has a 
long lifetime, bounded from below by about 10~ 2 sec. If the LOSP is charged or coloured it has 
a significant probability of being stopped in the detector, so that its decays can be observed. 
Studying the production and decay of the LOSP at the LHC could directly probe the freeze-in 
reaction that created the FIMP dark matter in the early universe. 
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